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We measure CP-violating asymmetries of neutral charmed mesons in the modes — > K~ K^ 
and D° — > 7r~7r"'" with the highest precision to date by using D° — » K~'k~^ decays to correct detector 
asymmetries. An analysis of 385.8 fb~^ of data collected with the BaBar detector yields values of 
ag|f = (0.00 ± 0.34 (stat.) ± 0.13 (syst.))% and a^p = (-0.24 ± 0.52 (stat.) ± 0.22 (syst.))%, which 
agree with Standard Model predictions. 

PACS numbers: 14.40.Lb, 13.25.Ft, 11.30.Er 



Charge-parity {CP) violation in decays of charmed mesons at levels as large as 1% has not yet been experi- 
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mentally ruled out [l| , and at this level would be evidence 
of unknown physical phenomena 0, H- The CP-even de- 
cays K^K+ and 7r~7r+ [3] are Cabibbo sup- 
pressed, with the two neutral charmed mesons, and 
I?", sharing the final states. CP-violating asymmetries 
in these modes are predicted to be C'(0.001%-0.01%) in 
the Standard Model of particle physics [5| , yet have not 
been measured precisely due to limited sample sizes and 
relatively large systematic effects Q. 

We search for CP violation in decays of charmed 
mesons produced from charm-quark pairs in the reaction 
e'^e^ — > cc by measuring the asymmetries in the partial 
decay widths, F, 



'■CP 



r{D" K-K+) - T{D" K+K- 
r(L)0 -> K-K+) + T{W 
T{D^ 7r-7r+) - r(^" - 



K+K- 



r(i:i" 



r(D" 
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In this construction, a^p, h = K,tt, includes all CP vio- 
lating contributions, direct and indirect 0]- The presence 
of direct CP violation in one or both modes would be sig- 
naled by a non-vanishing difference between the modes, 

Precise quantification of asymmetry in I?°-flavor as- 
signment, called tagging, has long been considered the 
primary experimental challenge in these measurements. 
We develop a new technique for measuring and correcting 
this asymmetry using only the recorded data. However, 
forward-backward (FB) asymmetry in cc production may 
be more significant at the center-of-mass energy of e"'"e~ 
collisions in BaBAR, « 10.6 GeV. This production 
asymmetry will create a difference in the numbers of re- 
constructed and D° events due to the FB detection 
asymmetries coming from the boost of the center-of-mass 
system (CMS) relative to the laboratory. 

The production asymmetry has two physical compo- 
nents. Interference in e"'"e~ — > cc as mediated by either 
a virtual 7 or a virtual Z° contributes at the percent 
level at this energy, and is well understood. In addi- 
tion, asymmetries induced by higher-order QED effects 
are expected to have polar angle dependence and to peak 
sharply in the forward and backward directions [?[. Al- 
though well-considered for /u-pair production , the pre- 
cise shape of this contribution for D production is not 
known. 

We use a data sample corresponding to an integrated 
luminosity of 385.8 fb~^ collected with the BaBAR detec- 
tor [i] at the PEP-II e+e" coUider at SLAC. The pro- 
duction vertices of charged particles are measured with 
a silicon-strip detector (SVT), and their momenta are 
measured by the SVT and a drift chamber (DCH) in 
a 1.5T magnetic field. Information from a Cherenkov- 
radiation detector, along with energy-deposition mea- 
surements from the SVT and DCH, provide K-t: discrim- 
ination. 



We analyze neutral D mesons produced from D*+ 
D^TT^; the charge of the tt^, a low momentum (soft) pion, 
indicates the flavor of the D° at production. To correct 
for asymmetry in this flavor tag, we measure the rela- 
tive detection efficiency for soft pions in recorded data 
using the decay D° ~* K^tt^ with (tagged) and without 
(non-tagged) soft-pion flavor tagging. The only detector 
asymmetry present in reconstruction of the signal modes 
is due to the tagging tt^, since the CP final states are 
reconstructed identically for and 

We reconstruct the four decay chains K^tt^; 
D*+ D°TT+, D° K-TI+- D*+ D^-Kt, £>" 
K-K+; and D*+ D°Tr+, D° 7r-7r+. We require 
D'^ candidates to have center-of-mass momenta greater 
than 2.4 GeV/c, which removes almost all B decays. Each 
daughter must satisfy a likelihood-based particle- 
identification selection and must have at least two po- 
sition measurements in each of the z and </> coordinates 
of the SVT. We require nf candidates to have a lab mo- 
mentum greater than 100 MeV/c and at least six position 
measurements in the SVT. 

For h ~ if , TT, we accept candidates with an invari- 
ant mass 1.79 < ruhh < 1.93GeV/c^ and, for final states 
with a TTg, an invariant mass difference 0.140 < Am < 
0.152 GeV/ c^, where Am = m/i^^^ — ruhh- For each D'^ 
candidate, we constrain the h'^h~ tracks to originate 
from a common vertex; for applicable final states, we also 
require the and tTs to originate from a common vertex 
within the e"'"e~ interaction region. We select candidates 
for which the probability of the vertex fit of the two 
Z?" daughters is greater than 0.005. For the KK and tttt 
modes, final asymmetries are calculated using events for 
which the polar angle of the 13° momentum in the CMS 
with respect to the beam axis satisfies | cos6'^^^| < 0.8. 

We statistically separate signal from background in 
the selected events by calculating signal weights based 
on an optimized likelihood function [10|. The likelihood 
function is composed of probability density functions 
(PDFs) that are fitted to the mass distributions using 
the maximum likelihood technique. For the non-tagged 
sample, a one-dimensional PDF is fitted to the rriKir 
distribution; for the tagged samples, two-dimensional 
PDFs are fitted to the rrihh and Am distributions. Two- 
dimensional PDFs are used for the tagged samples to ac- 
count for possible asymmetries in the background from 
correctly reconstructed decays with a misassociated 
TTs candidate; this background category peaks in irihh 
but does not peak in Am. The PDFs in this analysis are 
nearly identical to those used in an analysis of the de- 
cay if+TT" iJTj], since the signal shapes and back- 
ground sources are very similar. Although the PDFs are 
motivated by studies of simulated events, all of the shape 
parameters are varied in the fits to recorded data. Our 
selection of PDFs is treated as a source of systematic un- 
certainty. Because the signal shape is indistinguishable 
for D'-* and distributions, we use the same signal PDF 
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FIG. 1; Invariant mass distributions of the KK final state 
tagged as (a) and (b) D", and the tttt final state tagged as 
(c) D'^ and (d) D". Distributions of data (points with error 
bars) in the signal region 0.1434 < Am < 0.1474 GeV/c^ are 
overlaid with fitted PDFs (dashed line, shaded areas). The 
white regions under the central peaks represent signal events, 
the light gray misassociated vr^ events, and the dark gray 
remaining nonpeaking background. The data are shown over 
ranges extended beyond the fitted regions to illustrate the 
physical background shapes. 

to describe both flavors of a mode and fit it to them si- 
multaneously to reduce statistical uncertainties. The KK 
and TTTT invariant mass distributions for and Z?°, with 
fitted PDFs overlaid, are shown in Fig. [T] This analysis 
is sensitive only to ratios of D'^-signal yields to ZJ^-signal 
yields, and not to absolute yields, so the final results are 
relatively insensitive to the exact forms of the PDFs. 

The decay K~Tr~^ is chosen as a calibration 

mode because it provides an easily reconstructed inde- 
pendent sample with high statistics. However, detector 
asymmetries in reconstruction of the D'^ final state can- 
not be ignored (see Fig.[2^a,b)). These must be corrected 
to isolate the soft-pion asymmetry. 

Using the non-tagged Ktt sample, we produce a map 
of the relative reconstruction efficiency between and 
in this final state in terms of the momenta of both 
daughters, shown by components in Fig. [2lja,b). For 
each daughter, we consider the momentum magni- 
tude and polar angle in the lab with respect to the beam 
axis; these components are correlated. The daughters 
are, however, factorizable from one another. By consid- 
ering the normalized product of the K and tt efficiency- 
map components, wc obtain a four-dimensional relative- 
efficiency map for correcting D'^ —>■ K^n^ relative to 
if+TT". The presence of prompt decays not 
originating from a Z?*+ in the non-tagged sample extends 
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FIG. 2: Ktt efficiency-map components obtained from the 
non-tagged D^' daughters (a) K and (b) tt, and (c) tts ef- 
ficiency map with (d) statistical errors from the tagged Kn 
sample. Maps are produced from the ratios of candidate num- 
bers of D° to D°. 



the kinematic boundaries of the map but does not oth- 
erwise affect it. 

This Ktt map is used to weight the D° candidates in 
the slow-pion tagged Ktt sample, eliminating asymme- 
tries due to the D^/D^ daughters. Because all charm 
production is subject to the same production asymme- 
tries, these are simultaneously removed from the tagged 
Ktt sample by this correction. After the weights have 
been applied, the remaining asymmetry in the sample is 
due to the relative soft-pion efficiency. 

We produce a map of the relative soft-pion efficiency in 
terms of the pion-momentum magnitude and polar angle 
in the lab (Fig. [HJc)). Charm production is azimuthally 
uniform, and cf) is found to be uncorrelated with other 
momentum variables. Therefore, the 4> dependence is 
accounted for by an integrated scale factor. The uncer- 
tainties shown (Fig. Hfd)) are due to the statistical un- 
certainties in the sample yields. Signal-mode yields 
are weighted with this tTs map to correct for the soft-pion 
tagging asymmetry. The signal modes (with remaining 
production asymmetries) can thus be analyzed for evi- 
dence of CP violation. In Table HJ we list the raw and 
post-correction yields for the calibration and signal sam- 
ples in this analysis. In calculating these corrections, his- 
togram bins near kinematic boundaries with fewer than 
5,000 events are removed. 

CP violation would appear as an asymmetry in D'^ /D^ 
yields, independent of any kinematic variables. Because 
of the FB asymmetry in production, we calculate yield 
asymmetries as a function of cos6' = cos and de- 
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TABLE I: Signal yields in reconstructed modes. Listed uncertainties are statistical only. Corrections are applied only to D' 
samples, but all post-correction samples are restricted to the phase space of the correction map. 



Raw yields Post-correction yields 

Final state Corr. used D° W 

Kn 3,363,000 ± 6,000 3,368,000 ± 6,000 none — — 

Knns 705,100 ± 1,000 703,500 ± 1,000 if tt map 633,300 630,100 
KKtvs 65,730 ± 340 63,740 ± 330 tt, map 65, 210 63, 490 
TTTTTT, 32,210 ± 310 31,930 ± 310 tt^ map 31,900 31,760 



compose these into even and odd parts. We define 

n±(rn^ f)\ - ^ po {±\ COS e\) ~ n {±\ COS 9\) 

^ ' nDoi±\cose\)+n^oi±\cos9\) ^ ' 

acp = acp{cos6) « (a+(cos6') -fa" (cos 6*)) /2 (4) 

aFB{cos9) « (a+(cos6') - a-(cos6')) /2, (5) 

where n and rij^o are the numbers of signal events for 
and 1)° after applying the weights discussed above, 
acp is the even component and aFB(cos0) the odd com- 
ponent. Eqs. |4] and [5] are approximate as second-order 
terms in have been omitted. The even part, repre- 
senting CP-violating effects, would provide evidence of a 
uniform yield asymmetry. The odd part represents the 
production asymmetry, including higher-order QED con- 
tributions. From the several values of acp obtained as a 
function of | cos 6*1, we obtain a central value from a 
minimization. 

We consider three sources of systematic error to be 
significant. One source is the choice of PDFs used to 
describe the signal and background distributions, which 
affects the statistical background subtraction. We esti- 
mate this systematic uncertainty by substituting differ- 
ent background shapes in m and Am and an alternative 
two-dimensional signal shape in the fits to the tagged 
samples. Another source is the binning choices made 
and dependences in the TTs-efficiency correction. We esti- 
mate the size of this uncertainty by varying the number 
of bins and the required number of events per bin in his- 
tograms used to calculate efficiencies, and by adding a 
(j) dependence to the efficiency correction. We find the 
largest uncertainty here arises from the particular choice 
of binning in the TTj-efficiency map. Because the system- 
atic uncertainty in applying the TTs-efficiency correction 
is the same for both modes, we evaluate its size using 



TABLE II: Summary of systematic uncertainties. 
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FIG. 3: CP-violating asymmetries in (a) KK and (b) tttt, 
and forward-backward asymmetries in (c) KK and (d) tttv. 
In (a) and (b), the dashed lines represent the central values 
and the hatched regions the la intervals, obtained from 
minimizations. 

the larger signal sample. Finally, we consider the proce- 
dure for extracting acp- We vary the binning and the 
accepted range of | cos 611; the largest uncertainty comes 
from the latter. All other sources of systematic uncer- 
tainty are highly suppressed because the final states are 
reconstructed identically for and P)". We summa- 
rize the contributions to the total systematic uncertainty 
in Table |TT1 The smaller tttt sample size influences the 
calculation of its systematic uncertainty. 

For KK, we measure ag^ = (0.00 ± 0.34 (stat.) ± 
0.13 (syst.))%. For tttt, we measure aJJ. = (-0.24 ± 
0.52 (stat.) ± 0.22 (syst.))%. Statistical uncertainties of 
0.1% in the tTj correction have been included in the final 
statistical uncertainty values. The even and odd asym- 
metries for each mode as a function of | cos6'| are shown 
in Fig. [3l We conclude from the minimizations in 
Fig. [3l^a,b) that there is no evidence of CP violation in 
either of the Cabibbo-suppressed two-body modes of D° 
decay. This result is in agreement with Standard Model 
predictions. It also provides a new constraint on theories 
beyond the Standard Model W\, some of which predict 
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significant levels of CP violation in these modes. The 
asymmetries observed in Fig. |3l[c,d) represent the two 
Standard Model asymmetries discussed. Although an 
exact prediction of these forward-backward asymmetries 
does not exist, the observed values are consistent with 
expectations. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BaBAR. The collaborat- 
ing institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), CEA and CNRS-IN2P3 
(France), BMBF and DFG (Germany), INFN (Italy), 
FOM (The Netherlands), NFR (Norway), MIST (Rus- 
sia), MEG (Spain), and STFG (United Kingdom). Indi- 
viduals have received support from the Marie Gurie EIF 
(European Union) and the A. P. Sloan Foundation. 



* Deceased 

t Now at Tel Aviv University, Tel Aviv, 69978, Israel 

* Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

^ Also with Universita della Basilicata, Potenza, Italy 
^ Also with Universitat de Barcelona, Facultat de Fisica, 



Departament ECM, E-08028 Barcelona, Spain 
[1] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 
94, 122001 (2005); B. Aubert et al. {BABAR Collabora- 
tion), Phys. Rev. D 71, 091101 (2005). 
[2] Y. Grossman, A. L. Kagan and Y. Nir, Phys. Rev. D 75, 
036008 (2007). 

[3] G. Burdman and I. Shipsey, Ann. Rev. Nucl. Part. Sci. 
53, 431 (2003). 

[4] Unless otherwise indicated, particle types and decay pro- 
cesses imply also their charge conjugates. 

[5] F. Buccella, M. Lusignoh, G. Miele, A. Pugliese and 
P. SantoreUi, Phys. Rev. D 51, 3478 (1995); S. Bianco, 
F. L. Fabbri, D. Benson and I. Bigi, Riv. Nuovo Cim. 
26N7, 1 (2003). 

[6] W. M. Yao et al. (Particle Data Group), J. Phys. G 33, 
1 (2006). 

[7] Photon-box amplitudes interfering with the single- 
photon amplitude and Bremsstrahlung amplitudes in- 
terfering among themselves create a forward-backward 
asymmetry in e"'"e~ — » cc. 

[8] F. A. Berends, K. J. F. Gaemers and R. Gastmans, Nucl. 
Phys. B 63, 381 (1973); R. W. Brown, K. O. Mikaehan, 
V. K. Cung and E. A. Paschos, Phys. Lett. B 43, 403 
(1973). 

[9] B. Aubert et al. {BABAR Collaboration), Nucl. lustrum. 

Meth. A 479, 1 (2002). 
[10] M. Pivk and F. R. Le Diberder, Nucl. lustrum. Meth. A 

555, 356 (2005). 
[11] B. Aubert et al. {BABAR Collaboration), Phys. Rev. Lett. 

98, 211802 (2007). 



